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1.  Introduction 


Considered  is  a capital  budgeting  problem  with  linear  decision 
variables  that  can  be  either  continuous  or  integer  and  where  some  or  all 
of  the  associated  cash  flows  are  random  variables  that  may  be  statisti- 
cally dependent. 

It  is  formulated  as  a problem  of  linear  programming  under  risk 
for  which  mainly  three  approaches,  stochastic  programming,  linear  pro- 
gramming under  uncertainty  and  chance-constrained  programming,  have  been 
developed. ^ 

Stochastic  linear  programming,  an  approach  originally  suggested 

by  Tintner  [32]  is  primarily  concerned  with  the  probability  distribution 

2) 

of  the  objective  function's  value. 

Linear  programming  under  uncertainty  or,  as  it  is  sometimes  termed, 
stochastic  programming  with  recourse,  was  suggested  independently  by 
Beale  [ 4]  and  Dantzig  [10].  The  basic  innovation  of  this  approach  is  to 
amend  the  problem  to  allow  the  decision  maker  the  opportunity  to  make 
corrective  actions  after  the  random  variables  are  observed  and/or  pay 
penalties  for  constraint  violations. 

Charnes,  Cooper  and  Symonds  [ 8]  have  presented  a third  approach 

3) 

called  chance-constrained  programming  in  which  each  constraint  must  be 
satisfied  with  a certain  tolerance  probability,  i.e.,the  fulfillment  of 
the  constraints  is  not  required  to  be  guaranteed  as  in  linear  programming 
under  uncertainty. 

For  the  formulation  of  objectives  in  chance-constrained  program- 
ming approaches  different  criteria  have  been  discussed,  the  principal  ones 
outlined  in  Charnes  and  Cooper  [ 7]. 
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2.  A stochastic  capital  budgeting  problem 


In  this  paper,  the  capital  budgeting  problem  is  considered  to  be 
adequately  represented  by  a chance-constrained  programming  model  with  an 
objective  function  that  includes  the  expected  value  and  standard  deviation 
of  a firm's  horizon  value. 

The  application  of  chance-constrained  programming  to  capital  bud- 
geting problems  has  been  discussed  before.  See  for  example  Arzac  [ 1], 

Byrne,  Charnes,  Cooper  and  Kortanek  [ 3],  Byrne,  Cooper,  Charnes,  Davis 
and  Guilford  [ 5],  Hillier  [15],  Naslund  [24,25],  Naslund  and  Whinston  [27], 
Robertson  [28],  Struve  [31]  or  Spetzler  [30], 

The  main  purpose  of  this  paper  is  to  show  that  a stochastic  capital 
budgeting  problem  with  continuous  and/or  integer  decision  variables  can  be 
formulated  as  a convex  chance-constrained  programming  problem  that  can  be 
approximated  by  an  ordinary  (integer  or  noninteger)  linear  programming 
problem. 

Furthermore,  we  will  show  that  the  proposed  procedure  allows  the 
determination  of  an  optimal  cash  reserve  and,  in  addition,  the  considera- 
tion of  decision  opportunities  that  deal  with  an  unexpected  surplus  or 
deficit  in  periodic  capital  budgets  as  well  as  of  penalties  involved  if  the 
probability  constraints  do  happen  to  be  violated. 


2.1  Basic  formulation 

Let  us  consider  a capital  budgeting  problem  in  which  we  assume  that 

a firm  chooses  those  physical  and  financial  projects  which  maximize  its 

horizon  value  subject  to  physical  and  financial  constraints  in  distinguished 
4) 

periods.  Furthermore,  it  is  assumed  that  (a)  interdependences  between 
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projects  and  (b)  the  capital  market  conditions  can  be  expressed  in  an 
(integer  or  noninteger)  linear  programming  format^  such  that  the  capital 
budgeting  problem  can  be  formulated  mathematically  as: 

T 


maximize  z = E c.x 

(la) 

t 1=1 

subject  to  E g .x  = -D 

i=l  C1  1 c 

(t-l,..,T) 

(lb) 

x^  e X 

(1-1,.., T) 

(lc) 

0 < x1  < 1 

(1-1,.., T) 

(Id) 

where 

= vector  of  decision  variables,  each  representing  the  fraction  of  a 
project  (physical  capital  projects  as  well  as  borrowing  and  lending 
opportunities)  started  in  period  i, 

c^  = vector  representing  the  horizon  value  of  post  horizon  cash  flows 
associated  with  one  unit  of  x^, 

g j * vector  representing  the  net  cash  flow  in  period  t associated  with 
one  unit  of  x^, 

Dt  = the  total  amount  of  net  cash  flows  in  period  t associated  with  pro- 
jects realized  prior  to  the  start  of  the  planning  process, 

and  where  the  constraints  x^eX  (i=l,..,T)  can  be  expressed  in  an  (integer 
or  noninteger)  linear  programming  format. ^ 

The  elements  of  the  vectors  c_^  (i=l,..,T)  and  gfi  (t-1, . . ,T; i-1, . . ,t) 
and  the  (t-l,..,T)  are  assumed  to  be  either  constants  or  random  variables 
that  might  be  statistically  dependent.  If  some  elements  of  the  c^  or  even 
some  elements  of  the  or  some  of  the  Dt  are  random  variables,  the  hori- 
zon value  z is  also  a random  variable.  The  same  is  true  for  nt  (t-l,..,T), 


f 
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the  net  cash  flow  in  period  t,  defined  by 

t 

n.  = E g,.^,  + D (t-l,..,T)  (2) 

i=l 

if  some  elements  of  the  or  some  of  the  are  random  variables. 

In  this  case,  the  problem  formulation  (1)  must  be  regarded  as  an 
incomplete  formulation  of  the  decision  situation  as  it  does  not  reflect 
the  decision  maker's  attitude  towards  the  probability  distributions  of  the 
horizon  value  and  the  periodic  net  cash  flows.  It  may  be  completed  by  spe- 
cifying rules  that  allow  the  determination  of  solutions  for  the  decision 
problem  which  accord  with  the  decision  maker's  preferences. 

2.2  Decision  rules 

We  consider  here  only  rules  for  decision  situations  where  we  assume 
that  the  decision  variables  must  be  assigned  values  while  the  values  the 
random  variables  will  take  on  are  still  unknown  (zero-order  decision  rules). 
For  capital  budgeting  problems  attempts  to  allow  for  the  probability  distri- 
bution of  random  variables  have  been  concentrated  on 

a)  imposing  probability  constraints  on  periodic  net  cash  flows  while  maxi- 
mizing the  expected  value  of  the  random  objective  variable  (see  Hillier 
[15]  or  Naslund  [25])  or 

b)  introducing  an  objective  function  which  includes  not  only  the  expected 
value  of  the  random  objective  variable  but  a measure  for  the  dispersion 
of  its  probability  distribution.^ 

Both  approaches  have  been  extensively  discussed  by  their  authors 
with  regard  to  the  maximization  of  a decision  maker's  utility  function. 

As  for  the  choice  of  an  appropriate  dispersion  measure  in  the  second 
approach,  comparative  discussions  and  model  calculations  have  been  mainly 
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concentrated  on  the  variance,  semi-variance  and  standard  deviation  of 
the  objective  variable’s  probability  distribution  (see  Markowitz  [21], 

Bey  [ 2],  Hamada  [13],  Mao  [22,23],  Robertson  [28]  or  Van  Horne  [33,34]). 

In  this  paper  we  combine  both  approaches  as  we  assume  a decision 
maker  that  allows  for  the  probability  distribution  of  the  random  periodic 
net  cash  flows  nt  (t=l,..,T)  and  the  random  horizon  value  z in  the  sto- 
chastic capital  budgeting  problem  (1)  by  preferring  decisions  which  are 
obtained  according  to  the  following  rules: 

1)  In  each  period  t (t=l,..,T)  the  equality  constraint  (lb)  should  hold 
for  the  expected  value  E(nt)  of  the  periodic  net  cash  flow  while  the 
probability  that  a periodic  deficit  exceeds  a specified  amount 
should  be  lower  than  a specified  risk  level  (l-at),  i.e.  , 


E(nt)  = 0 


(3) 


Prob  {n  > -b  } > a 
t - t - t 

2)  A possible  realization  z'  of  the  horizon  value  z,  defined  by 

z'  = E(z)  - B J V(z)  with  8^0  (4) 

with  the  expected  value  E(z)  and  the  standard  deviation  / V(z)  of  the 
probability  distribution  of  z,  should  be  as  high  as  possible.  For  a 
discussion  and  determination  of  6 in  a probability  context  see>for 
example, Hi llier  and  Lieberman  [16]. 

While  the  second  rule  focuses  on  the  standard  deviation  of  the 
horizon  value's  probability  distribution  as  a dispersion  measure  to  be 
used  in  the  objective  function  of  the  programming  model,  it  is  easy  to 
show  that  other  measures  like  the  variance  or  semi-variance  of  the  distri- 


bution (as  preferred  by  Markowitz  [21])  can  be  used  within  the  frame- 
work of  the  proposed  procedure  as  well. 


2.3  A complete  formulation 


By  introducing  the  decision  rules  into  the  programming  formula- 
tion (1)  of  the  stochastic  capital  budgeting  problem,  it  can  be  formu- 
lated as  a chance-constrained  programming  problem: 


T J T 

maximize  z'  = E(  Z c.x.)  - B/  V(  Z c.x. ) (5a) 

i=l  i=l 

t 

subject  to  E(  Z g .x.)  = EC-0,.)  (5b) 

i=l 

t 

Prob  {(  Z g^Xj+D,.)  t -bt*  - at  (t=l,..,T)  (5c) 

xi  e X (i=l,..,T)  (5d) 

0 >_  XjL  >_  1 (i=l,..,T)  (5e) 


The  first  step  in  solving  this  problem  is  to  reduce  it  to  a deter- 
ministic equivalent  form,  i.e,,to  determine  deterministic  equivalent  forms 
of  the  probability  cjnstraints. 

3.  A deterministic  equivalent  form 

Denote  by  st  a "situation"  k in  period  t,  representing  a possible 

Ic  1c  lc 

simultaneous  realization  [g^t , • • ,g^t .D^ ] of  the  random  elements  of  g^t  for 

i=l,..,T  and  and  by  sT+1  a situation  representing  a simultaneous 
k k 

realization  [c^,..,^]  of  the  random  elements  of  ct  (t=l,..,T).  Let  St  be 

the  set  of  all  possible  situations  in  a period  t (t*l,..,T)  or  the  horizon 

/ 

period  T+l,  respectively. 

Furthermore,  denote  by  s^  a sequence  j of  situations  over  the 
planning  period  such  that 

sJ  = [s-j,sj,..,s^+1]  with  s^  e St 


(t*l,..,T+l) 
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and  let  S be  the  set  of  all  possible  sequences. 

For  the  moment,  it  is  assumed  that  S includes  only  a finite  small 
number  q of  possible  sequences  and  that  these  sequences  and  their  proba- 
bilities pJ  (j=l,..,q)  are  known.  The  possible  net  cash  flows  in  the 
various  periods  and  the  possible  horizon  values  are  then  determined  by 
(a)  the  value  of  the  decision  vectors  xfc  (t=l,..,T)  and  (b)  the  possible 
realizations  of  situation  sequences. 


3.1  The  probability  constraints 


Let  us  now  consider  the  possible  net  cash  flows  in  a single  period  t 

for  a specified  value  of  the  decision  vector  x^  (i=l,..,t).  Denote  them  by 

n^  (j=l,..q)  and  regard  them  as  the  possible  realizations  of  a random 

variable  n . Furthermore,  denote  by  E(n  ) the  net  cash  flow  in  period  t 
tx  ’ tx 

of  an  "expected  situation  sequence"  representing  the  expected  values  E(gti), 
E(Dt)  and  E(c^)  of  g^,  f>t  and  c^  with  i=l,..,T. 


Define 


-j  fE(ntx>  - nL 
tX  ( 0 

4 Kx  ~ E(ntx) 
tX  ( 0 


nJtx  i E(ntx)  | 
otherwise  * 


n^  > E(n  ) ) 
tx  tx  I 

otherwise  j 


Let  us  assume  that  the  functional  form  of  the  probability  distri- 
bution of  n is  known  and  that  the  fractiles  of  this  distribution  are 
tx 

8) 

completely  determined  by  its  mean  and  variance.  Let  F(u)  denote  the 
cumulative  distribution  function  of  the  standardized  variable 
u - (n  -E(n  )//  V(n  ) . Define  u by  the  relationship  F(u  )«a  . 

L a L.  A La  UL  v»L  L 


Then  the  deterministic  equivalent  form  of  the  probability  constraint  in 

9) 

period  t for  a specified  decision  vector  can  be  formulated  as 


E(ntx)  ‘ Uat/  V(ntx)  - -bt 


However,  if  the  functional  form  of  the  probability  distribution 
of  nfcx  is  not  known,  Tchebychev's  inequality  (or  the  extension  of 
Camp-Meidel  if  ntx  can  be  regarded  as  unimodal,  see  Duncan  [11])  still 
yields  an  upper  bound  for  ua(.  This  bound  is  a very  conservative 

measure.  However,  it  will  be  shown  that  the  model  calculations  yield 
some  information  about  the  functional  form  of  the  probability  distribu- 
tion of  the  net  cash  flow  associated  with  the  optimal  decision  vector. 
This  information  might  be  used  (a)  to  adjust  the  assumption  about  the 
functional  form  or  (b)  to  calculate  the  "true"  a^,  i.e.  the  which 
refers  to  the  used  uat  and  the  updated  information  about  the  functional 
form  of  the  net  cash  flow's  probability  distribution. 

Whatever  the  probability  distribution  of  n happens  to  be,  as  its 
expected  value  E(ntx>  has  been  restricted  to  zero  in  (5),  (7)  can  be 
rewritten  as 


£ -J«^)2  + <<x>2>  i (V“«)2 


(8)  defines  a convex  upper  bound  for  the  absolute  deviation  of  the 
possible  n^x  (j»l,..q)  from  E(ntx).  It  restricts  the  set  of  feasible 
values  of  the  decision  vector  x^  (i=l,..,t)  in  order  to  satisfy  the 
probability  constraints  in  (5). 
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As  the  . llue  of  the  decision  vector  is  not  known  explicitly,  the 
values  of  the  k^  (j=l,..,q)  and  k^  (j»l,..,q)  in  (8)  are  not  either. 
However,  they  can  be  determined  simultaneously  with  the  determination  of 
an  optimal  decision  vector  by  introducing  variables  k (j=l,,.,q)  and 
k^  (j=l,..,q)  with  the  possible  realizations 

kt^  : ktx  V xi  e x (i=l,,.,t) 

vi  (9) 

k?  ! ktx  V Xi  E X (i=l,..,t) 

into  the  programming  model  (5).  The  deterministic  equivalent  form  of  the 
probability  constraint  in  period  t for  an  unspecified  decision  vector  can 
then  be  formulated  as 

£ Pj((k"j)2  + (k^j)2)<  (bt/uat)2  (10) 

j=l 

From  an  economic  point  of  view,  these  variables  could  be  regarded 
as  single  period  borrowing  and  lending  opportunities  for  unexpected  cash 
flows  in  period  t.  In  this  case,  borrowing  rates  rt  and  lending  rates  r^ 
might  apply  such  that  the  borrowing  or  lending  of  a cash  unit  in  period  t 
would  affect  the  net  cash  flow  in  period  t+1  by  Rt=(l+rt)  or  R*=(l+r*), 
respectively.  While  we  will  follow  this  interpretation  when  formulating 
the  deterministic  equivalent  form  of  the  chance-constrained  programming 
problem,  the  variables  could  be  considered  alternatively  as  representing 
penalties  involved  if  the  constraints  (5b)  do  happen  to  be  violated. 
However,  in  this  case  the  penalties  might  have  to  be  paid  in  period  t 
rather  than  in  period  t+1,  i.e.,the  model  formulation  would  have  to  be 
changed  accordingly. 


i 
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3.2  Objective  function 


■ 


A similar  procedure  can  be  used  to  reduce  the  objective  function 

of  (5)  to  a more  tractable  form.  Let  us  define  a random  variable  w bv 

x } 


the  relationship 


T T 

w = E c x - E(  £ c.x*) 
i=l  i=l  1 1 


where  x^  (i=l,..,T)  represents  a specified  value  of  the  decision 
vector  x^  (i=l,..,T).  From  (11)  it  follows  that  w^  has  an  expected  value 


E(w  ) — 0 and  a variance 
x 


V (w  ) = V(  £ c.x  ) 
i=l  1 1 


Now,  let  us  denote  the  possible  realizations  of  w^  by 

T * T • * • * T * , 

f E(  £ c.x  ) - £ c^x.  if  c-?x.  < E(  £ c.x.) 

-j  1 i=l  1 1 i=l  1 1 1 1 ~ i=l  1 1 


j=l» • • »q 


otherwise 


T . T T 

(£  c^x*  - E(  £ c x*)  if  cj.x*  > E(  £ c.x*)') 

“ 1=1  “ 

0 otherwise  J 

where  j marks  the  relevant  situation  sequence  s^  (j=l,..,q).  Then  we  get 

V(w  ) = £ pJ((w"j)2  + (w+j)2)  (14) 

J-l  X 

By  proceeding  as  described  above  for  the  determination  of  the  variance 
of  a periodic  net  cash  flow  we  introduce  variables  w ^ (j=l,..,q)  and 
w ^ (j“l,..,q)  with  the  possible  realizations 


I 


into  the  programming  model.  Let 


y2  = V(w  ^+w+^)  and  v = E(  Z c x.) 

i=l  :L 

Then  we  can  replace  the  objective  function  of  the  chance-constrained 
programming  formulation  by  the  equivalent 


maximize  z ' = v - |Jy 
T 

subject  to  E(  Z c x ) - v = 0 (16) 

i=l  1 1 

Z pj((w_j)2+(w+j)2  - y2  < 0 

j=l 

where  the  quadratic  constraint  defines  a convex  bound  for  the  set  of 
feasible  solutions  of  the  decision  problem. 

3. 3 Model 

The  incorporation  of  the  variables  k ^ , w ^ and  w+^  with 
t=l,..,T  and  j=l,..,q  into  the  programming  model  requires  the  explicit 
consideration  of  each  possible  situation  in  period  t (t=l,..,T)  and  the 
"horizon  period"  T+l,  i.e.  the  consideration  of  all  q situation  sequences 
over  the  planning  period.  (This  procedure  is  similar  to  the  one  used  in 
stochastic  programming  with  simple  recourse  [38]). 

The  deterministic  equivalent  form  of  the  chance-constrained  pro- 
gramming problem  can  then  be  formulated  mathematically  as: 


Maximize  z'  = v - 6y 


subject  to 


= -E(D1) 


E(gll)xl 

8nxi  + kiJ  - k?  = 


4 

^ (p  j (k~  j ) 2+pj  (k*j  ) 2)  < (bx /u^) : 


+ J’v'pjRt-ikt-i+p3Rt-ik^i> 

.,E18tixi  - "t -l\-l  + Rt-lkt-x  + ktj  - ktJ 


-E(Dt)  (t 


0 

t' 


z (pJ(k~j)2+pj(k+j)2)  < (b  /u  )2  (t= 

i=l  c c - t at 


_^E(c1)xi  + - v 

T 

- R^k,^  + R^k^  - v + w j - w+j 


E c~!  x 


i=l 


ii 


0 (j* 


E (p^(w  ^)2+p-^  (w+^)2) 


j=l 


2 

y <0 


x±  e X (i-1 

0<xt  <1  (i-1 

kt'5,ktJ,v,w"J,w+J,y  -°  (J*1 
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1 

As  the  nonlinear  constraints  in  (17)  are  all  based  on  quadratic 
separable  functions,  they  can  easily  be  approximated  by  a grid  lineari- 
zation.  In  this  case,  efficient  (integer  or  noninteger)  linear  pro- 
gramming algorithms  can  be  used  to  solve  the  problem. 

It  is  easy  to  show  that  the  deterministic  equivalent  form  of  the 
chance-constrained  programming  problem  can  be  adapted  to  more  complex 
situations,  i. e. , situat ions  where  the  decision  maker's  opportunity  set 
for  dealing  with  an  unexpected  periodic  surplus  or  deficit  consists  of 
several  distinguished  single-period  or  multi-period  borrowing  and 
lending  opportunities.  Furthermore,  a technically  similar  formulation 
would  allow  the  consideration  of  penalties  that  have  to  be  paid  if  a 
deficit  in  a period  t exceeds  specified  amounts  or  if  the  probability 
constraints  do  happen  to  be  violated. 

As  an  example,  consider  a situation  where  several  distinguished 
single-period  borrowing  and  lending  opportunities  are  available.  Denote 
by  y ^ and  vectors  that  represent  the  borrowing  and  lending  opportu- 
nities in  a period  t within  the  situation  sequence  j.  Assume  borrowing 
and  lending  rates  to  be  represented  by  vectors  and  R*  . Furthermore, 
assume  that  the  capital  market  conditions  for  these  opportunities  can 

be  expressed  in  a (integer  or  noninteger)  linear  programming  format  re- 

-1  +1 

presented  by  the  formulation  y ,y  Je  A . 

Then  the  constraints  in  (17)  that  refer  to  a period  t=2,..,T 
could  be  replaced  by  the  following  formulation: 
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+ Y 


-J 


R+  Y+j  1 

Rt-iVi} 

* -E(Dt) 

R+  +j  + k-j  _ k+j 

t-lYt-l  t Kt 

= -H5 
t 

(J-i,...q) 

5 (pJ(k"j)2+pj(ktj)’?> 
j-1 

1 (bt/uat)2 

- - k;j  + 

= 0 

(j=l, . . ,q) 

y-i  Y+J 
Yt  .Yt 

e A 

0=1,.., q) 


4 . 2 Consideration  of  a variable  cash  reserve 

The  bt  (t=l,..,T)  in  the  constraints  (5c)  can  be  regarded  as 
restrictions  on  the  amount  of  cash  that  is  easily  available  from  exter- 
nal sources  for  the  compensation  of  "unexpected"  deficits  in  the  perio- 
dic net  cash  flows,  i.e.,  their  negative  deviation  from  the  expected 
values.  The  formulation  (5)  does  not  allow  the  consideration  of  inter- 
nal resources  for  this  purpose.  This  is  due  to  the  assumption  usually 
made  in  chance-constrained  programming  approaches  for  capital  budgeting 
problems  (see »for  example, Naslund  [25])  that  the  decision  variables 
represent  "contracts"  that  cannot  be  changed  during  a period  for  com- 
pensating unexpected  deficits. 

However,  the  possibility  of  determining  an  optimal  periodic 
cash  reserve  from  internal  resources  can  be  incorporated  into  the 
programming  formulations  (5)  and  (17). 


(18) 
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Let  us  assume  that  external  and  Internal  resources  are  available 

for  compensating  unexpected  deficits  in  a period  t.  Denote  by  y a deci- 

t 

sion  variable  that  represents  the  amount  reserved  from  internal  resources. 
Then  (5b)  and  (5c)  in  (5)  can  be  replaced  by 


E<.E  - h ■ E<-Dt> 

i=l 


Prob  {(^gt^+Dj.)  > -(bt+yt)}>  ot 


yt  > ° 


(t=l,..,T) 

(t-l,..,T) 

(t-l,..,T) 


(19) 


For  reducing  the  probability  constraints  in  (19>  to  a deterministic 

* 

equivalent  form  we  define  variables  y by  the  relationship 

(t=l,..,T)  (20) 


= b + v 
t t ' t 


and  proceed  as  described  in  the  preceding  sections.  We  then  get 


uat  ^(pj(ktJ)2+(k*j)2)-  (yt)2<0  (t=l,..,T) 


(21) 


yt  - yt  >-bt  (t=i,..,T) 


The  quadratic  constraints  that  define  convex  bounds  on  the  set  of 
feasible  solutions  of  the  problem  are  based  on  quadratic  separable  func- 
tions, i.e. , their  introduction  into  the  programming  model  (17)  would  not 
affect  the  possibility  of  its  approximation  by  a linear  programming 
formulation. 


5.  Random  sampling  of  situation  sequences 


Suppose  now  that  the  number  q of  possible  situation  sequences  is 
large  and  that  it  is  not  feasible  or  economical  to  consider  all  sequences 
explicitly  in  the  model  formulation.  In  this  case,  the  decision  of  the 
decision  maker  might  be  based  on  a random  sample  of  sequences. 

Assuming  that  the  future  probability  distributions  of  the  random 
variables  and  their  interdependences  are  known,  a random  sample  of  situa- 
tion sequences  could  be  generated  by  a simulation  process.  However,  the 
available  data  about  future  events  is  often  restricted  to  a number  of 
sequences  usually  derived  from  the  past,  i.e.,from  time  series, etc.  Then 
this  set  of  sequences  is  usually  regarded  to  represent  a random  sample. 

However  a random  sample  of  situation  sequences  happens  to  have 
been  formulated,  its  use  affects  the  formulation  of  the  probability  con- 
straints of  the  problem.  If  the  probability  distributions  of  the  periodic 
net  cash  flows  (or  the  horizon  value)  can  be  regarded  as  normal  or  appro- 
ximately normal,  the  necessary  changes  in  the  formulation  of  the  deter- 
ministic equivalent  forms  of  the  probability  constraints  can  be  calculated 
in  a straightforward  manner. 

Let  us  consider  a probability  constraint  in  the  form  (7)  and 
neglect  for  simplification  the  indices  x and  t such  that 

E(n)  - u^  ✓ V(n)  > -b  (22) 

_ 2 

From  the  sample  mean  m and  the  sample  variance  s we  may  compute 
(m-u^s)  as  an  estimate  of  (E(n)-u  1/  V(n)  ).  In  repeated  sampling  from  a 
stable  distribution  this  estimate  will  be  approximately  normally  distri- 
buted (see  Hald  [12])  about  (E(n)-u  / V(n))  so  that  we  cannot  be 


sure 
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that  (m-u^s)  as  computed  from  a single  sample  will  be  larger  or  equal 


(-b)  even  if  E(n)-ua/  V(n)  <►  -b  . To  be  "reasonably"  sure,  i.e.,vith  a 
specified  probability  p,  that  the  probability  constraints  will  hold,  the 
multiplication  factor  u in  (m-u  s)  has  to  be  replaced  by  a factor  u 

a a r j g 

such  that 

Prob  {(m-u  s)  > (E(n)-u  / V(n)>  > p (23) 

P — a — 


It  has  been  shown  elsewhere  (see  Hald  [12])  chat  the  normal  distri- 
bution that  approximates  the  probability  distribution  of  X=(m-u„s)  has 

D 


mean  M(A)=(E(n)-u0/  V(n))  and  variance 
p 


V(X)  = V(n)  ( i - 2I  ) 


(24) 


where  d denotes  the  sample  size  and  f*(d-l)  the  degree  of  freedom.  Now, 
let  F(u)  denote  the  cumulative  distribution  function  of  the  standardized 


variable  u=(X-M(X)) //  V(X)).  Define  u by  the  relationship  F(u  )=p  . 

P p 

Then  by  proceeding  in  the  usual  way,  the  probability  expression  (23)  can 
be  replaced  by 

M(X)-up/  V(X)  > (E(n)-uo/v(n)  (25) 


or 


./l  uj» 

-U.  - U 1/  -T  — -TTe  >~U 

P p V d 2f  — a 


(26) 


As  u , u , d and  f are  assumed  to  be  known,  11  can  be  calculated  from  (26) . 
pa  p 

If  the  functional  form  of  the  probability  distributions  of  the  net 
cash  flows  nfc  (t=l,..,T)  are  not  normal,  the  sampling  distributions  of 
means  and  variances  will  not  be  the  same  as  if  they  were  normal,  but  the 
computations  will  not  be  very  seriously  affected  unless  the  departure  from 
normality  is  very  marked  and  the  samples  are  small  (see  Duncan  [11]). 


If  the  functional  form  is  not  known  or  we  are  unwilling  to  make 
assumptions  regarding  it,  then  it  is  still  possible  to  derive  the  mean 
and  variance  of  X=(m-UgS)  from  the  sample  (see  Duncan  [11]).  But  no 
information  is  available  about  the  functional  form  of  the  probability 
distribution  of  X because  no  assumptions  can  be  made  about  the  functio- 
nal form  of  the  distribution  of  s.  In  this  case  one  might  apply  Tcheby- 
chev's  inequality  to  get  some  rough  limits  on  the  probability  variation. 

6.  Conclusion 

In  the  preceding  sections  we  were  considering  a capital  budgeting 
problem  with  continuous  and/or  integer  decision  variables  where  some  or 
all  elements  of  the  data  were  statistically  dependent  or  independent  ran- 
dom variables.  It  was  formulated  as  a chance-constrained  programming 
problem  that  allowed  the  explicit  consideration  of 

a)  decision  opportunities  dealing  with  a deficit  or  surplus  in  periodic 
net  cash  flows  such  as  the  accumulation  of  an  optimal  cash  reserve 
or  appropriate  borrowing  and  lending  opportunities  and  of 

b)  penalties  that  have  to  be  paid  if  periodic  deficits  do  occur  or  the 
probability  constraints  do  happen  to  be  violated. 

It  could  be  shown  that  the  chance-constrained  programming  problem 
can  be  approximated  by  an  ordinary  (integer  or  noninteger)  linear  pro- 
gramming problem.  The  proposed  procedure  is  based  on  the  explicit  consi- 
deration of  possible  realizations  of  the  random  variables  in  the  program- 
ming formulation  as  discussed  in  linear  programming  under  uncertainty  for 
stochastic  programs  with  simple  recourse.  The  resulting  increase  in  the 
size  of  the  programming  formulation  could  be  limited  by  basing  the  pro- 


■ 


U 


— 


, I |g  | 


cedure  only  on  those  possible  realizations  of  the  random  variables  that 
were  connected  with  a random  sample  of  situation  sequences,  i.e.,of 
possible  subsequent  periodic  realizations  of  the  random  variables  over 
the  planning  period. 
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Footnotes 


1.  For  a survey  see  Naslund  [25],  Vajda  [36]  or  Kail  [18]. 

2.  See  Sengupta  [29]  for  a comprehensive  discussion. 

3.  See  also  Charnes  and  Cooper  [ 6]. 

4.  See  Charnes,  Cooper  and  Miller  [ 9]  for  the  original  presentation  of 
a horizon  value  objective  function  and  Weingartner  [37]  for  a compa- 
rative discussion  of  models  using  the  horizon  value  or  the  present 
value  of  a firm  in  their  objective  function,  respectively. 

5.  See  Weingartner  [37]  for  some  specific  formulations. 

6.  Note  that  the  constraints  (lb)  are  formulated  as  equalities.  This  is 
due  to  the  assumption  that  the  set  of  bending  and  borrowing  opportu- 
nities in  each  period  t includes  all  available  alternatives  for 
cash  use. 

7.  See  Markowitz  [21]  for  the  original  presentation  though  in  the  context 
of  a portfolio  investment  problem. 

8.  See  Hilller  [15]  for  a discussion  of  relations  between  the  probability 
distribution  of  the  random  variables  g^  (i=l,..,t)  and  Dt  and  the 
probability  distribution  of  their  linear  combination  n and  the  con- 
ditions under  which  the  probability  distribution  of  n can  be  regar- 
ded as  normal  or  at  least  (by  some  version  of  the  Central  Limit 
Theorem)  approximately  normal. 

9.  See  Hillier  and  Lieberman  [16]. 

10.  See  Hillier  [14]. 
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